Inherent Structures in models for fragile and strong glass 



Luca Leuzzi and Theo M. Nieuwenhuizen 
Universiteit van Amsterdam 
Valckenierstraat 65, 1018 XE Amsterdam, The Netherlands 
(printout: February 1, 2008) 

O 

' An analysis of the dynamics is performed of exactly solvable models for fragile and strong 

, glasses, exploiting the partitioning of the free energy landscape in inherent structures. The results 

are compared with the exact solution of the dynamics, by employing the formulation of an effective 
' temperature used in literature. Also a new formulation is introduced, based upon general statistical 

considerations, that performs better. Though the considered models are conceptually simple there 
is no limit in which the inherent structure approach is exact. 
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1. INTRODUCTION 



o 

1^ I The characteristics of a glassy system ||l|,|| arises from the complex topography of the multidimensional function 
representing the collective potential energy that gives rise to a non-trivial partition function and thermodynamic 
potential. In this picture, at low enough temperature where vibrations are minimal, the spatial atomic patterns in 
, crystals and in amorphous systems share the common basic attribute that both represent minima in the potential 
^ ' energy function describing the interactions. The presence of distinct processes acting on two different time-scales 
means that the deep and wide local minima at and below the glass transition temperature Tg are geometrically 
1 organized to create a two scale length potential energy pattern. Tg depends on the cooling procedure and it is usually 
^ [ determined as the temperature at which the viscosity of the glass former reaches the value of 10^'^ Poise. 

In the present paper we investigate, using the inherent structure approach, an exactly solvable model glass that 
shows all the basic features of real glasses |^. The model is built by processes evolving on two different, well separated 
time scales, representing respectively the a and (3 processes taking place in real glassy materials. The slow a processes 
^ ' represents the escape from one deep minimum within a large scale valley to another valley. The fast (3 processes, 
instead, are related to elementary relaxations between neighbouring minima inside the same valley. We consider here 
all kinds of /3 processes as equivalent, since the characteristic time-scales on which they are evolving are in any case 
much shorter than the time-scale of the a processes (i.e. the observation time). 

In the general case, decreasing the temperature, the free energy local minima can, in principle, be split into smaller 
' local minima. But if we can assume that they maintain their identity in spite of this splitting, we can set a one-to-one 
T-H , correspondence between local minima and inherent structures i. e. between the minima of the free energy and 

' the ones of the potential energy. Actually, such a splitting is not even present in the two dynamical models presented 
' here after, making the correspondence clearer. 
!^ In this paper we will see to which extent such a scheme, widely used in numerical simulations P-p|,p|-p^ , applies to 

-««.^ our analytically solvable model. We will compare it with the exact dynamic solution, achieved without any partitioning 
of the configuration space. 

In section ^ we do introduce the two kinetic models and we give the description of their statics and of their Monte 
Carlo dynamics. In section |^ we develop the inherent structure approach for the dynamics of such models and we 
'■^ define two different inherent structure effective temperatures mapping the dynamics into a thermodynamic frame (in 
^ [ section 3.3); one definition follows the literature about numerical simulations P p^, the other exploits the analytic 



THE MODELS AND THEIR PROPERTIES 
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' 1. Hamiltonian and its constraint 

We present two dynamical models, having the same statics, but different dynamics bringing to the behavior of a 
fragile glass in one instance and to the behavior of a strong glass in the other one. The version describing a system 
relaxing like a fragile glass was introduced in |l^ and widely studied in 
Both models are described by the following local Hamiltonian: 
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where N is the size of the system and {xi} and {Si} are continuous variables, the last satisfying a spherical constraint: 
Sf = N. We call them respectively harmonic oscillators and spherical spins. K is the Hooke elastic constant, H is 
an external field acting on the harmonic oscillators, J is the coupling constant between {xi} and {Si} on the same site 
i and L is the external field acting on the spherical spins. A separation of time scales is introduced by hand: the spins 
represent the fast modes and the harmonic oscillators the slow ones. We assume that the {Si} relax to equilibrium 
on a time scale much shorter than the one of the harmonic oscillators. From the point of view of the motion of the 
{xi}, the spins are just a noise. To describe the long time regime of the {xi}, in ||^, we did average over this noise 
by performing the computation of the {Si} partition function, obtaining an effective Hamiltonian depending only on 
the {xi}, that determines the dynamics of these variables. Using the saddle point approximation for large N we got: 
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We define, then, the effective Hamiltonian 'Hcs{{xi}) 
ration of {xi}. We find 
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This can also be written in terms of the internal energy U{{xi}) and of the entropy Sep{{xi}) of the equilibrium 
processes (i.e. the spins): 
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The function U is actually the Hamiltonian averaged over the spins and Sep is the entropy of the spins. 

In we studied the model characterized by a constraint on the phase space, introduced for the fragile glass case to 
avoid the existence of the single global minimum, and implementing a large degeneracy of the allowable lowest states. 
The constraint is taken on the {xi}, thus concerning the long time regime. It reads: 
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where mg is a model parameter. It is a fixed but arbitrary, strictly positive constant. This constraint applied to the 
harmonic oscillators dynamics is a way to reproduce the behavior of good glass formers. We imposed a Monte Carlo 
dynamics [ p^ , p^ satisfying this constraint and coupling the otherwise non-interacting {2;^} in a dynamic way. As we 
saw in the system exhibits a Vogel-Fulcher-Tammann- Hesse (VFTH) relaxation ||2^, characterizing a fragile glass. 

To model a strong glass, instead, we will also consider the same model Hamiltonian but without imposing any 
constraint and making use of a different Monte Carlo dynamics. We will show later in this paper that this dynamics 
displays an Arrhenius relaxation near zero temperature. In this case we have a strong glass, as it happens for similar 
models, e. g. the oscillators model and the spherical spins model where exactly the same dynamics is applied. 
The new point of the present model is that now both fast and slow processes occur. 

To shorten the notation we define the modified "spring constant" K and "external field" H: 
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K = K - 
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H = H 
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We stress that K and H are actually functions of the {xi] themselves (through mi and m2 that occur in w). We also 
define the constant 



D = HJ + KL. 
Using the definitions ( 2. 10] ) it is useful to note that 

HJ + KL = HJ + KL = D. 



(2.11) 
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2. Statics at heat-bath temperature T 

The partition function of the whole system at equilibrium is: 



Z{T)= j VxVSe^^[-m{{xr},{S.})]5\Y, 
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The new object that appears in the exponent is the configurational entropy 
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It comes from the Jacobian expX of the transformation of variables Vx dmidm2, (see ( |2.3| )). We can compute the 
large A^-limit of this partition function using once again the sadd le po int approximation. The saddle point equations 
are found minimizing the expression between square brackets in (^.13|) with respect to mi and m2. This yields 
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The form of the solutions rni{T), m2{T) is quite complicated because each of these equations is actually a fourth 
order equation, but they can be explicitly computed. In terms of the equilibrium values Wik we find the following 
expression for the equilibrium free energy: 
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This is the statics both for the model with the constraint ( |2.9D , as long as the temperature exceeds the Kauzmann 
temperature, and for the o ne w ithout it. Indeed for the fragile glass case at T <Tq, when the constraint is reached, 
the saddle point equation ( ^.1(S| ) becomes 7772 — TOi = toq, no matter what the temperature of the thermal bath is. In 
this work, however, we will limit ourselves for the fragile glass to cases where T is slightly larger then Tq, and for the 
strong glass case to temperature slightly above zero. 



3. Dynamics 



The dynamics we apply to the system is a parallel Monte Carlo dynamics firstly introduced in ||lj|. The thus 
obtained dynamical model composed by the simple local Hamiltonian ( p.l] ) and such a dynamics has the benefit of 
being analytically solvable. 
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In a Monte Carlo step a random updating of the variables is performed {xi x'^ = Xi + Vi / \/ N) where the {r^} have 
a Gaussian distribution with zero mean and variance A^. We define x = Ti-dx'^) — TCdxi}) as the energy difference 
between the new and the old state. If x > the move is accepted with a probability W{l3x) = exp(— else it 
is always accepted {W{/3x) = 1). The updating is made in parallel. It is the parallel nature of the updating that 
allows the collective behavior leading to exponentially divergent time scales in models with no interactions between 
particles such us ours. A sequential updating would not produce any glassy effect. This dynamics may induce glassy 
behavior in situations where ordinary Glauber dynamics ]l7t would not. In our model the parallel dynamics mimics 
the presence of interactions between atoms in realistic glasses, where a high internal cooperativeness is present. For 
different examples of dynamics implying non trivial collective behavior the reader can look, for instance, at the n-spin 
facilitated kinetic Ising model [ p^Jl9| or at the kinetic lattice-gas model |2^j2l|] . 

In a Monte Carlo step the quantities Nmi = Xi and Nm2 — J2i ^'''^ updated. We denote their change by yi 
and y2, respectively. Their distribution function is, for given values of mi and 7712, 

P{yi,y2\mi,m2) = y]J-^^==e-''''/(2A^) (5 (^a;', -^x, -yij -5 ( ^ x' ■ - ^ a;,^ - j 



Neglecting the variations of mi and m2 of order A'^ /N we can express the energy difference as Q 

x^^y2-Hyi, (2.20) 

In terms of x and y — yi the distribution function can be formally written as the product of two other Gaussian 
distributions: 

Piyi,y2\mi,m2)dyidy2 ^ dx p{x\mi,m2) dy p{y\x,mi,m2) 

dx ( {x-xf \ dy ( {y-y{x)f \ 

^ exp — exp — (2.21) 



^/IIK^ V 2A^ y ^2^ V 2A, 



where 



■^A^k/2, A^ = A^K^{m2-ml) + A^K^(mi~H/ky, (2.22) 

nii—H/K x — x A^(m2 — m^) 

m2 — rrii + imi ~ H / K \ m2 — mf + Imi — H / K j 



y{x) = rn,-H/K ^ ^x - x ^ ^ AHm2 - mj)^ ^ ^ ^^.23) 



a. Dynamics of the fragile glass model 

To represent a fragile glass the dynamics that we apply to the system is a generalization of the analytic treatment 
of Monte Carlo dynamics introduced in . As noted in , also in this generalized case the dynamical model with 
a contrived dynamics can be analytically solved. As we saw in in the long-time domain, the dynamics looks quite 
reasonable with regard to what one might expect of any glassy system and the system exhibits a VFTH relaxation. 
We repeat here the main steps of the implementation of this dynamics (for a more extended presentation see j^]). 

We let A^, the variance of the random updating {r^}, depend on the distance from the constraint, i.e. on the whole 
{xi} configuration before the Monte Carlo update: 

A^it) ^ 8[m2(i) - m?(t)] [ ^ V (2.24) 

\m2(t) - mf(t) - mo J 

where B, ttiq and 7 are constants. In particular 7 is an exponent larger than zero that appears in the VFTH- 
like relaxation law of the model, when T decreases towards some critical temperature Tq (identified in |3| with the 
Kauzmann temperature): 
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where A{ is a constant depending on the system's parameters. In other models p^ p^ , p3t the variance was kept 
constant. The same will also be done for the strong glass version of the model in the next subsection. 

For what concerns the exponent 7 we saw in |^ that it generates different dynamic regimes for 7 > 1,7 = 1 and 
< 7 < 1; the situation 7 = 1 remains model dependent even in the long time limit. We will stay in the following in 
the regime for 7 > 1. 

The nearer the system goes to the constraint (i.e. the smaller the value of m2 — — toq); the larger the variance 
becomes, implying almost always a refusal of the proposed updating. In this way, in the neighborhood of the 
constraint, the dynamics is very slow and goes on through very seldom but very large moves, that can be interpreted 
as activated processes. When the constraint is reached the variance A^ becomes infinite and the system dynamics gets 
stuck. The system does not evolve anymore towards equilibrium but it is blocked in one single ergodic component 
of the configuration space. At large enough temperatures, the combination m2{t) — — tuq will remain strictly 

positive. The highest temperature, Tq, at which it can vanish for t — > 00, is identified with the Kauzmann temperature 

In p] the dynamics was expressed in terms of two combinations of mi and m2. The first one, defined as 

fii = ^-mi. (2.26) 
K 

represents the distance from the instantaneous equilibrium state. By instantaneous equilibrium state we mean that 
H and K depend on the values of mi and TO2 at a given time t. For t —^ 00, at the true equilibrium, one has fii = 0. 
The second dynamical variable is defined as the distance from the constraint ( ^.91 ): 

fi2 = m2- ml - mo- (2.27) 

When fi2 = the constraint is reached. This will happen if the temperature is low enough (T < Tq) and the time 
large enough. Tq is the highest temperature at which the constraint is asymptotically {t —^ 00) reached by the system. 
Above Tq ordinary equilibrium will be achieved without reaching the constraint. The temperature is, then, too high 



for the system to notice that there is a constraint at all on the configurations, and this implies (see (2.16)) 



where 



lim ^i2it) = MT) = - TOO > , (2.28) 



lim K {mi{t,m2{t);T) = K^{T) ^ K {mi{T),m2{T)) . (2.29) 

t^OO 



Below To the system goes to configurations that become arbitrarily close to the constraint, and then stay there 
arbitrarily long. Note that, by definition of Tq we can write 

Too = — (2.30) 

KM) 

Solving the equations of motions, for fixed parameters (aging setup), we find, to the leading orders of approximation 
for large times, the following behavior for /i2 

H2it) ^ Tj- (2.31) 

[log(tAo) + clog(log(tAo))]'/'' 

where c = 1/2 since in this paper we only look at the regime for T > Tq. The co nstant depends on the parameters 



of the model and on the temperature; it is of order one. The solution ( 2.31 ) is valid in the aging regime, where 



to ^ i ^ Teq{T). Indeed, when t ~ Teq{T) ~ exp (A/(T — To)) the "distance" ^2 becomes 

B 



M2 



[( 
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T-To) 





(xT-Tq , (2.32) 



as it should be. 

We also introduce another variable that will be useful later on, namely the difference between fJ.2(t) and its asymp- 
totic, equilibrium, value 'p2(T): 
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(2.33) 



where, using (2.3C), ^jl^{T) comes from ( p.28 ), valid, in the fragile case, when T > Tq. When t — ^ cx) is, by definition, 
S^i2 = 0. 

The dynamical behavior of /ii depends not only on the temperature (above or below Tq) but also on 7 being greater, 
equal to or lesser than one. With respect to the relative weight of /zi and /i2 we can identify different regimes, where 
the solution has different behaviors What is of our interest here is the regime of T > Tq and 7 > 1, where 
fJ-i{t) <^ ^J■2it) and a unique effective thermodynamic parameter can be properly defined in various independent ways 



b. Dynamics of the strong glass model 

We now analyze the simple case without constraint on the configuration space and where A^, the variance of the 
randomly chosen updating {r,} of the slow variables {xi}, is a constant. This dynamical model can also be seen as 
the limit for mo and 7 ^ 1 of the preceding one. We also mention that the case with J = L = is the model of 
harmonic oscillators studied in [p^ , p3| . 

In the fragile glass case we studied a different version of such a dynamics for two particular combinations of the 
variables mi an d m? . Here we will keep the same notation. The first variable is thus defined, starting from the saddl e 
point equation ( 2.15 ), as the deviation from the instantaneous equilibrium state and is formally equivalent to ( 2.26 ). 

The second variable is defined as 

fi2 = iTi2~ ml. (2.34) 



When T = from equation (2.16) we know that /i2 = 0. Indeed at T = the system reaches its minimum 

X, = V* . (2.35) 

For simplicity we limit ourselves to a choice of the interaction parameters such that D — HJ + KL > and K > 0, 
for which this is the global minimum. In appendix A we derive the equations of motion for ni and H2 and we solve 
them for temperature equal to and slightly above zero and long times, in the aging regime. In this time regime fj,i 
comes out to be much smaller than ^2- /^i oc The solution for ^2 comes out to be, at the leading order 

,.it) . (2.36) 



The difference between /i2(^) and its asymptotic value is now: 



6^.2{t) ^ ,2{t)-ll2{T) ^ - ^ (2.37) 



where /i2(T') comes from ( 2.16| ) and 

KM) - \m^K{m,{t),m2{t);T) = + ^j^^ + ^ d{D + J^)^ + ^ ^'-^^^ 
For t 00, (5/i2(^) 0. 

At low temperature, the relaxation time for the slow processes depends on the temperature following an Arrhenius 
law: 

req(r) « exp (^^ j ; (2.39) 

A..^. ,2^40, 
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4. Two temperature thermodynamics 



Before going on we recall here that we are able to introduce effective parameters in order to rephrase the dynamics 
of the system out of equilibrium into a thermodynamic description (for a review see |p3|| ). 

In 1^ we got through different methods the following expression for the effective temperature in the regime for 
T > To as a function of the interaction parameters of the model and of the time evolution of its observables: 

Te{t)^ K{mi{t),m2{t))[mo + fi2{t)] ■ (2.41) 

Since we will use one of these methods in the next section to map the IS dynamics into an effective thermodynamic 



parameter we shortly recall this particular derivation of (2.41). Knowing the solution of the dynamics at a given time 



t a quasi-static approach can be followed by computing the partition function Z,. of all the macroscopically equivalent 
states at the time t. In order to generalize the equilibrium thermodynamics we assume an effective temperature 
and an effective field H^, a nd substitute the Boltzmann-Gibbs equilibrium measure by exp(— 7ieff({a:i}, T, iJe)/Te), 



where TicS is given in (2J) and the true external field H in it has been substituted by the effective field He- As we 
get the expression of the "thermodynamic" potential Fe = —T^ log Zp as a function of macroscopic variables mi^2 
and effective parameters, we can determine Tg and Hf, minimizing F^, with respect to mi and m2 and evaluating the 
resulting analytic expressions at toi^2 = ^1,2(0- 

The partition function of the macroscopically equivalent states is: 



Ze{mi,m2]Te,He) = j Vx exp 



-7^Hoff({xJ,r,i/e) 



SiNmi-y^x,) S{Nm2-22xl)^ (2-42) 



From this we build the effective thermodynamic potential as a function of Tf, and He, besides of T and H, where 
the effective parameters depend on time through the time dependent values of mi and m2, solutions of the dynamics. 
They are actually a way of describing the evolution in time of the system out of equilibrium. The free energy 
Fe = —TelogZe is minimized with respect to mi and m2. Then their time dependent values are inserted, yielding 



Fe{t) = U (mi(i), m2(t)) - TSep {mi{t), m2{t)) - Te{t)I (mi(i), m2{t)) + [H - H e{t)]N mi{t) , (2.43) 



with 



Teit) = k (mi{t),m2{t))[mo + ^l2{t)] , 

He{t) ^ H ~ K imi{t),m2it)) fii{t) . (2.44) 

where Tg is the so called effective temperature. U is the internal energy of the whole system, Sep is the entropy of the 
fast or equilibrium processes (the spherical spins) while T is the entropy of the slow, " configurational" , processes (the 
harmonic oscillators). The last term of Fe replaces the —HNuii occurring in U (see eq. (p??|)) by —HeNmi. U , Sep 
and X are 'state' functions, in the sense that they depend on the state described by T, Te, H and, if needed. He- In 
the framework where only one relevant effective parameter Te stays, these functions do not depend on the path along 

which its value has been reached. 

As we saw in |^ for the VFTH relaxing model at T > Tq and in appendix A for the Arrhenius relaxing case ( A.25| ), 



( A.35| ), the effective temperature alone is enough for a complete thermodynamic description of the dominant physical 



phenomena {He — H). The introduction of He becomes important only for second order corrections in 6122- 



3. INHERENT STRUCTURE APPROACH 



The characteristics of a glassy system can be represented by means of a multidimensional potential energy function 
with a complex topography. The spatial patterns of atoms in crystals and in amorphous systems, at low temperature, 
represent, then, minima in th e po tential energy function describing the interactions 



In the case of the model (2.1) all the complex chemical properties of real glass formers are not present, but 
nevertheless the system exhibits several aspects of their complex features (described in Q), indicating that our simple 
model is complicated enough for what concerns the description and the comprehension of the basic long time properties 
of a glass. 

In a real glass the presence of distinct processes (acting on different time-scales) can be obtained from a careful 
analysis of the relaxation response function above Tg. We limit ourselves to a two time-scale approach. This means 
that the deep and wide local minima at and below Tg are geometrically organized to create a two scale length potential 
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energy pattern. As a consequence the system shows a and [3 processes. The a processes represents the escape from 
one deep minimum within a large scale valley to another valley. This escape requires a lengthy directed sequence of 
elementary transitions producing a very large activation energy. Moreover, the high lying minima between any two 
valleys, among which the system is making a transition, are very degenerate. This implies a large activation entropy 
for the inter-basin transition. f3 processes are instead related to elementary relaxations between neighboring minima 
{intra-basin dynamics). 

We stress that in our models we put together all kind of (3 processes in our short time-scale, since they are in any 
case much shorter than the observation time considered. 



1. Decomposition of the partition function: introduction of inherent structures 

In this point of view an approximate approach to the problem is to divide the complicated multidimensional 
landscape of the (potential) energy in structures formed by large deep basins and to describe the dynamics of the 
processes taking places as intra-basin and inter-basin . 

More precisely one can define an Inherent Structure (IS) as that basin behind an actual configuration of the system 
evolving in time at some temperature T that is the minimum of the potential energy reached in an instantaneous 
quenching by the method of steepest descent. 

The introduction of IS's allows, at low enough temperature (T < Tg), a decomposition of the partition function 
into an IS part, connected to the zero temperature landscape corresponding to the configurations of the system at 
temperature T, and a part connected to the thermal excitation of the configurations in a single minimum. 

The probability that an equilibrium configuration at T — belongs to a basin associated with an IS structure 
with an energy density in the interval [e, e + de] is 

V{e,T)de oc exp (-/?7V [e - rsc(e) -t- /„(e, T)]) de (3.1) 

where Sc(e) is proportional to the logarithm of the number of IS's existing at the energy level e and /«(e, T) is the free 
energy of the configurations inside an IS at energy e (related to a temperature T system) . To derive the distribution 



(3.1) in this form the approximation is made that /„ is computed as the average over all the IS's of energy e. This 
means that, by assumption, the shape of a basin depends only on its energy level and on the temperature. Enough 
below Tg the further approximation can be made, that fy{e,T) ~ fy{T), because fiuctuations inside one IS are small 
The shape of the basin depends then only on the temperature. All the internal (vibrational) states of any IS 
have the same (vibrational) free energy fv at given T. We anticipate, however, that in the present study we will not 
carry out such an approximation for our models. 

IS dynamics is significant, i.e. significantly represents the actual dynamics of the system at finite T, provided that 
there is a one to one correspondence between IS's and real minima at finite temperature and provided that these IS's 
are visited with the same frequency with which the corresponding finite T minima are visited. 



2. Inherent Structure approach in the Harmonic Oscillator-Spherical Spin model 



As we will see, the model (2.1) is built in such a way that every {xi} configuration is an inherent structure. 
Indeed, at a given {xi} configuration at finite T, the {Si} are fast variables and they contribute to the energy and 
to the other observables as a noise depending on temperature. If we take away this contribution we do not actually 
change the configurations of the minima of the slow variables. In the case of the system without constraint on the 
configuration space, nor contrived dynamics, any {xi} configuration is an inherent structure. For what concerns the 
constrained model, instead, certain configurations are not allowed. Moreover the presence of the constraint ( |2.9| ) 
produces (entropic) barriers higher than in the other case to get from a certain IS to a different one. That just means 
that the dynamics through the inherent structures is even slower in the fragile glass case that in the strong glass case. 

First of all wc have to define the steepest descent procedure for the model. We start performing the minimization 

of 

N 

H + Xj^Sf-^N, (3.2) 



where Ti. is the Hamiltonian (2.1) of the model and where wc implemented the spherical constraint '^^Sf — N hy 
using the Lagrange multiplier A. 
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To get rid of the contribution of the spins, i.e. to get rid of the fast modes, we minimize ( p.2[ ) with respect to the 
{S,}. We get 



S.. 



(min) _ PjJXj + L) 

2X 



Inserting this values for Si and solving for A we, then, find 

A = -Wis 



where mi and m2 are defined in (2.3) and 



= \/ Pm-i + 2JL'mi + 



Using (3.4) the minimum {Si\ configuration for a given set of {xi\ is, thus, given by 

-,(min) JXi L 



s 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



Finally the expression (3^) becomes 



N 



K 

yTO2 - Hmi - Wis 



(3.7) 



that is the energy function of the inherent structures. Consequently the partition sum over inherent structures is 
defined by 



Zis = j X'x exp [— /J'His] = j dmidm2 exp {X — PHis) 



(3.8) 



Due to the minimization any explicit dependence on T disappears with respect to the effective Hamiltonian ( |2.5| ). In 
( ^.2D , integrating over th e sp ins, instead of minimizing with respect to them, we also had an entropic term for the fast 
processes {Sep given in ( |2.S| )) and a slightly different internal energy {Nw instead of Nwie). Carrying out steepest 
descent the entropic term vanishes (only the minimal configuration is taken into account) and the inherent structure 
energy has no explicit dependence on the temperature. The configurational entropy for IS's comes from the Jacobian 
of the transformation of variables Vx = e^dmi dm2: (see (2.3) and ( 2.14 )). It is the same of the finite T case, since 
any allowed configuration {xi\ is also an IS. 
The static average of TLis is given by 



E^(T) = n,Mt:l{T)) 



(3.9) 



where Tn^i^.i'^) ^'^s the solutions of the saddle point equations that we get in the IS case to compute (3.8), in the limit 
of large N. The equations are 



D 



1 JiK-jywis) J Ki 

■(is)N2 



_(is) 



(m\ ') = —{Jm\ ' +L) = ^ 



where we define 



JIL 

Wis 



(3.10) 
(3.11) 

(3.12) 



with Wis from (3_^). The combination HisJ + KigL — HJ + KL = D is, again, simple, as in ( 2.12 ). 
In the case at finite T the static partition function (p.l3|) was 



Z = j dmidm2 exp (X — /37icff), 



(3.13) 
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with TicS defined in (2.5) and T in (2.14). The two saddle point equations are different from ( ^.ISj ) and ( |2.16| ) valid 



in the realistic case, giving thus different results: m\ 2 7^ Tni^2- We note explicitly that m\ 2 depends on T even in 
the IS case. 

Comparing the expressions so far obtained with those appearing in the exponent of the probability distribution 



(3.1) we identify the configurational entropy Nsc with X, as defined in (2.14), and the rest with 



iV(e + U) - {{x^}) = His ({a:.}) + F„ {{x,}) (3.14) 

where, as already told, His {{xi}) is the IS internal energy and from the difference HcS {{xi}) — "Wis {{xi}) = Fy{{xi)} 
the thermal free energy of one IS turns out to be: 

F,^^\og(^^^^^)-N{w-w..) . (3.15) 



where w is defined in (2.4) and in (3.5). Notice that it explicitly depends on the parameters mi and m2 of the 



IS, whereas in literature it is often assumed to be a constant (harmonic approximation |4y5|,|8| |l l| . 

3. Effective temperature in the IS's approach 

a. Expansion of the dynamical energy 

A possible way of defining an effective temperature, sometimes used in literature, for instance in the study of 
Lennard- Jones interacting spheres (^H] and in the study of the random orthogonal model ||l^, is to compare the 
time dependent out of equilibrium mean internal energy with the equilibrium mean internal energy expression at a 
temperature Te 7^ T. The out of equilibrium mean internal energy is built taking the dynamics of a system out of 
equilibrium at temperature T and repeating it many times starting from different initial conditions. A statistical 
ensemble of trajectories is constructed in this way. At any given time t the configurations that each sample is visiting 
are found. The energy Tiis averaged over the ensemble of different trajectories is, 

E^J'\t)= {n^s)t = N^m^it) - NHmi{t) - N^/j'^m2{t) + 2JLmi{t) + (3.16) 



^Y"^2'' - iViJmf ^ - N\J ,Pmf'^ + 2 JLmf ' + L"^ + A^i^,,(mf \ m^''^)J^2(0 + C{mf\mP)5^l2{tf (3.17) 
E^^\T) + NkUr4'\rnf'^)5m{t) + C{mf\m!i'^)5Mtf (3.18) 



where Sfj,2{t) = f^2{t) — M2(^) given by (2.37) in the Arrhenius case and by ( ^.33 ) in the constrained case for T >To 



The equilibrium IS energy i?cq ' (T) will be a different function of the temperature in the two dynamic versions of the 
model. The second order expansion will be needed only for the strong glass case and the expression for the factor 
C{nii'^\ 7712^'') is, in that case: 

= ■ (3-19) 

We can then take a system in equilibrium at a temperature Te, such that the configurations visited by the system 
at equilibrium are the same as those out of equilibrium at temperature T. This we call effective temperature. In other 
words, fixing i, is defined as the temperature at which the system at equilibrium would visit the same configurations 
visited by the system out of equilibrium at temperature T, with the same frequency. 

b. The effective temperature employed in numerical approaches: the fragile case 

Following the approach found in literature for numerical simulations, we can define a Tg}^'' through the 

matching of the equilibrium and the out of equilibrium IS internal energy: it is the one such that 

E^(T^^{t))^El^^^{t) (3.20) 
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For our model it is possible to work out an analytic expression for such a T^'l^\t), at least near the Kauzmann 
transition for the fragile glass case, linearizing in T — Tq. 

What we get is a parameter different from the thermodynamic effective temperature ( 2.41 ) that we got from three 
different approaches (including the Fluctuation-Dissipation Ratio) in jsj. 

For the fragile glass case we arc not able to derive any simple expression, of the IS energy ( 3.16 ), but we can in 



any case solve it exactly. The results are shown in figures 3.1-3.2 for a given choice of the values of the interaction 
parameters of the model. As one can see Tj:'^\t) comes out to be different from Te{t) at any time decade. 

As a matter of fact what we are comparing now with the average E^{t) is a function E^glT^'^'') of the effective 
temperature alone, while we know that out of equilibrium any proper thermodynamic function cannot simply depend 
on just one temperature as the thermodynamic function of equilibrium systems do p^ ]. It is not surprising, thus, 
that the two functions do not coincide. 



c. The effective temperature employed in numerical approaches: the strong case 



For the strong glass case it is possible to work out a simple analytic expression for the dynamical and the equilibrium 
IS energy. To do it we will expand near zero temperature up to second order in T. 

We underline that the thermodynamic effective temperature given in (2.41) is also the expression of the effective 
temperature for the system without constraint. What changes in that case is the time behavior of m2 — ml = /i2, 
that is now given by ( |2.36 ), and its limit at equilibrium (see ( 2.16| )). In this case, where the analytic treatment is by 
far easier, we can give a short explicit expression for T^'^\t): 



.(is) 



T 



KD 
D + .P 



2(D + J2) 



Here above terms of 0{T^) and 0{6fj,2it)^) cancel. This T^[^\t,T) is obtained from (|3.20D with 



E'iit)^ {H + Jf ^ T J^K 
N 2K 2 8D{D + J^y 



KD 



TTSfl2{t) 



(3.21) 



2(D + J2) ^ SD{D + J^f '^^^ ' %[D^,PY 
O (T^) + O {T^6fi2{t)) + O {T6fi2itf) + O {6fji2{tf) (3.22) 



If we expand (2.41) in the same way we get: 

KD 



T,= T + KSfi2{t) = T 



JK 

d1~p 



DJ^K^ 
2{D + .P)* 



Sfl2{tf 



-O (T^) + 0{TH^l2it)) + O {T5ti2{tf) + O {5ii2{tf) 



(3.23) 



As we see from the formulae above and from figures |3.3| and |3.4| , for a given choice of the parameter values, in the 
case with Arrhenius relaxation Tg and T^i'^ are very similar. Their difference is one order of magnitude less than in 
the model with contrived dynamics. 



d. A more fundamental definition of the IS effective temperature 



Here we propose an alternative way to identify an effective temperature that maps the dynamics between inherent 
structures into a thermodynamic quantity. We follow a quasi static approach using a partition sum, just as we did in 
the finite T case. The aim is to be able to define an effective thermodynamic parameter for the IS dynamics and to 
compare it with the given in (2.41). Following exactly the same approach we used in |3j (see section 2. 4), including 



the substitution of the real external field H with the effective one Hj,2' , we compute the partition function counting 
all the macroscopically equivalent IS's, through which the system is evolving in this symbolic dynamics, at a given 
time t. 
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10'" 10"" 10™ 10*" 10^'" 10°" 10'" 10'^" 10"" 10'"" 
FIG. 3.1. Effective tempera ture s vs. t at the heat bath temperature T = 4.00248, equal to the Kauzmann temperature. The 
constants in the Hamiltonian (2.1) are set to the foUowing values: K = J— 1, H — L = 0.1. The constraint constant is mo = 5. 



The upper curve shows the effective temperature got by matching out of equilibrium and equilibrium IS internal ener gy. T he 
one in the middle is the behavior of (2.41), for systems at finite T, and the lowest one is the IS effective temperature (3.29). 





lO-* 10'" 10" 10"" 10"-' lO-"" 

FIG. 3.2. The same effective temperatures, for the same choice of parameters as before are plotted for a different heat bath 
temperatureiT = 4.1. 



K 



(is) 



5{ Nmi ~^xA 5 Nm2 -^x"^^ 



Tj(is) — 



T. 



(is) 



^— log (m2 — m 



exp • 



[-P^^ [His ( 



TYis ( mi , m2 ; T, h''^^'' ) " t';;^2''I (mi , m2 ) 



(3.24) 

(3.25) 
(3.26) 



f^^e2 ~ ^/'^e'i' ^^'^ ^e2' parameters describing the behavior of the system going only through IS's. Minimizing 



.(is) 



r(is) 



the free energy 



(is) 



(is) (is) 

-T^2 logZe with respect to mi^2 we get: 

.(is) _ 



^e2 = ^is (mi, m2) [m2 - mf\ , 

r(is; 
'e2 



(3.27) 
(3.28) 
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By in serting the time dependent values of mi and m2 we now look at the time evolution of the effective temperature 
( |3.27 ) for lar ge ti mes, in the aging regime, and we compare it with the behavior of the thermodynamic effective 
temperature (2.41). 

For the dynamically constrained model, for t oo, T^2'' — > T (if T > Tq). When to ^ i < oo, h owever, the way 
the effective temperature approach the heat bath temperature is different from the behavior ( ^.41 ) of Tg, found in 
the case at finite temperature. For a comparison their first order expansions are: 



(is) 



e2 



T 



Te~T 



1 + T 



1 + T 



2D{l + Qt^D) 



20(1 + Q^D) 



(3.29) 
(3.30) 



where 



Xis.^(r) = lim k,,{mi{t),m2{t)-T) 

t—>-oo 

k^{T) = lim if (mi(t),m2(t);r) , 

t — *CXD 

n(-) = lini 



t^oo K-w- 

IS IS 



lim Q = lim — 

t^OQ t~^oo ^^11) ^yj 



-T/2Y 



(3.31) 



(3.32) 



(3.33) 



The time dependent variable 5ii2{t) is introduced in ( ^.33 ) in both cases (apart from the parameter influencing 
only the short times) while the coefficients in front of it are different at any temperature, including Tq. In the 

fragile case, thus, this second IS effective temperature does not coincide with T^-^ and it is much nearer, at any 

time, to (3.30). However, even if this t'^2^ is conceptually more properly chosen than the one defined matching out 



of equilibrium energy at temperature T and equilibrium energy at temperature Ti*'*'' , we still do not get the same 
parameter describing the finite T dynamics in a thermodynamic frame. The inherent structure approach gives thus 
a good approximation but is nevertheless never analytically correct. 

To show how good this approximation is we can take as an instance a certain realization of the model with given 



values of the "fields" and "coupling constants". We plot in figures |3.1| and |3J the behavior of T^'^\t), T^^' {t) and 
Tf,{t) at heat-bath temperatures equal to and just above the Kauzmann temperature. 

For the strong glass case we also expand for temperatures near to zero and for long time and we get: 



.(is) / 



T^^\t)^T + k,,5pL2{t) = T- 



KD 



5ll2{t) 



l5n2{t) 



-5ii2{tf 



+0{T^) + 0{T^Sii2it)) + 0{T5^2{tf) + O {5^2{tf) 
Te{t) - 2p + j2)3 ^^/^2W + ^ (^A^2(0') + 0{T^) + 0{T'Sf,2{t)) + 0{T5ii2m 

T^i\t) + ^^^^M^t)^ + ^(^') + o{TH^l2{1^) + o{T5^l2{tr) + o {5^,2{tf) 



(3.34) 

(3.35) 
(3.36) 



where 5^,2 is given by (2.37). 

The effective temperature mapping the dynamics of the system evolving at finite temperature T have the same 
behavior of Te2 in approaching the heat-bath temperature up to order T5^2{i) where they start deviating one from the 
other. For a quenching to zero temperature the two effective temperatures coincide. Moreover, due to the simplicity 



of the model, the IS effective temperature is equal to with T)^-^' given in (3.21) up to order 5fi2{t) hi time and 
up to order in temperature. 



(is) 
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10 10** 10" 10° 10'" 10" 10" 

FIG. 3.3. Time evolution of the effective temperatures at the heat bath temperature T — 0.0005 in the model with Arrhenius 
relaxation. The constants in the Ham iltonian (2.1) are set to the following values: K = J=1, H = L = 0.1. The lower curve 
shows the effective temperature (3.21) got by matching o ut of equilibrium and equilibrium IS internal energy. To order 5^2 it 

Second order differences are too small to appear in the plot. 

!.4l' 




coincides analytically with the IS effective temperature (3.35) 
The upper curve is the behavior of 



for systems at finite T. 



0.01 






10^ 10^ lO"* 10^ lO" lO' 

FIG. 3.4. The same effective temperatures, for the same choice of parameters as before are plotted for a different heat bath 
temperature: T = 0.001. Comparing the time scales of the two plots, we can clearly observe the decreasing of the Arrhenius 
relaxation time to equilibrium Teq that takes places rising the temperature. 



4. CONCLUSIONS 



In this paper we consider a model that owns all the basic properties of a glass, built by processes evolving on two 
well separated time scales, representing the a and (3 processes taking place in real glassy materials Q. The decoupling 
of time scales is fundamental for a generalization of equilibrium ther mod ynamics to systems far from equilibrium. 

We take into account two different versions of the model given by (^]^) . One leading to the description of a fragile 
glass having a non-zero Kauzmann temperature and the other one representing a strong glass. 

Using a particular Monte Carlo dynamics and developing it analytically, thus having the opportunity of probing it 
in more detail with respect to a numeric study, we found equations of motion that are in all respect those typical of 
glass relaxation. 

In the strong glass case we apply exactly the same parallel Monte Carlo dynamics used in [p^-^ 



23 1, finding an 
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Arrhenius relation between the relaxation time of the slow processes {xi} and the temperature. 

In the fragile glass case the model is provided with a constraint applied to the harmonic oscillator dynamics, i. e. 
to the slow processes dynamics, in order to reproduce the behavior of a good fragile glass former. In j^], by means 
of a Monte Carlo constrained dynamics, we identified the Kauzmann temperature with the one, Tq, at which the 
constraint is reached, asymptotically, for the first time in a cooling experiment from high temperature. There we 
showed how the thermodynamic phase transition |^5|] is characterized, that takes place due to the breaking of the 
ergodicity in the landscape of our model, rich of degenerate minima. 

In this work we carried out the inherent structure approach. In both dynamical models, decreasing the temperature, 
the free energy local minima do not split into smaller local minima, just like in the p-spin model in zero magnetic field 
|p6| , because every allowed configuration of harmonic oscillators is and stays an inherent structure at any temperature. 
Consequently we can set a one-to-one correspondence between the minima of the free energy and the ones of the 
potential energy (i. e. the inherent structures). Because of this exact correspondence the dynamics through inherent 
structures should be a valid symbolic dynamics for the real system, i.e. at a finite heat-bath temperature T. At least, 
it would significantly represents the actual dynamics if the inherent structures are visited with the same frequency 
with which the corresponding free energy minima at finite T are visited. 



We defined the steepest descent procedure for the model, that is the minimization of the effective Hamiltonian (2.5) 
appearing in the partition function with respect to the spins, i. e. the fast relaxing variables. Due to the minimization 
any explicit dependence on T disappears. The configurational entropy for inherent structures was computed from the 
logarithm of the Jacobian of the transformation of variables Vx —^ dmidm2, and thus it was the same of the exact 
finite T approach. In our models, then, any configuration of harmonic oscillators {xi} (for the fragile glass model 
every configuration allowed by the constraint), is also an inherent structure. Although the models we considered 
are conceptually very simple and without interactions, as compared to another approach proposed for systems with 
interacting discrete spins where the IS scheme breaks down p7t , our setup seems to be more physical since it is 
intimately based on time scale separation between fast and slow processes. A direct consequence of this time scale 
separation is that we encounter a both mathematically and physically well defined configurational entropy, whereas 
this observable suffers from principle difficulties in |2^]. 

(is) 

We can take a system in equilibrium at an effective temperature Te , such that the configurations visited by 
the system at equilibrium are the same as those out of equilibrium at temperature T. First we defined an effective 
temperature through the matching of the equilibrium and the out of equilibrium internal energy of the inherent 
structures (the one such that EiT^it)) = E^(t)). For the strong glass model this effective temperature almost 
coincides with the finite temperature Tg provided that the temperature at which the system is quenched is not too 
high (as far as terms of 0{TS^2{t)) are negligible they are equal). On the contrary, when the constraint is set and 

(is) 

the contrived Monte Carlo dynamics is applied, we found that the thus derived effective temperature Tg\ is quite 
different from the effective temperature that we were able to identify in the finite T dynamics. Therefore we proposed 
a new definition following a quasi static approach. In this way we computed the partition function counting all the 
macroscopically equivalent inherent structures, through which the system is evolving in this symbolic dynamics, at 
a given time t. Even though the result we get is much more similar to the finite T dynamics effective temperature 
(numerically speaking the difference is one order of magnitude smaller), yet it is analytically different, indicating 
that the inherent structure scheme can only be an approximation to what happens in the realistic dynamics of the 
system. As a consequence, also the derivation of out of equilibrium thermodynamic quantities (e.g. the configurational 
entropy) obtained making use of this approach could suffer of a systematic deviation from the exact result. 
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APPENDIX: STRONG GLASS DYNAMICS 

In this appendix we present the Monte Carlo dynamics of the observables /ii and /i2, functions of the slow relaxing 
harmonic oscillators {xi} through 
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in the case where the model (2.1) is not subjec ted to any constraint on its {xi} configurations 



Let us recall the definitions, given in section 2. 3 



H , 

fii = — — mi (A. 2) 
K 

fj,2 = m2- m\ (A. 3) 

In this notation the average and the variance (2.23) of the gaussian distribution p{x\mi,m2) (see ( ^.21 )) of the 
possible changes in energy during the Monte Carlo dynamics become: 



(A.4) 



We remember here that x is the difference (2.20) between the energy of the configuration proposed for the exchange 
and the energy of the actual configuration. A is fixed. 

To shorten the following expressions we also define the parameter 



A2 



V2A; V8(M2 + M!) 



(A.5) 



A{t) = J dx W{Px) p(x|mi,TO2) 
and the rate of change of the energy of the system: 

m^JdxxWmpixlm^m,) 

Defining the auxiliary function: 

/3^A, ' ^ 



where 



f{t) = a;/3exp —/3a; + 



erfc(a) 



rfc 



dz e 



A, 

— p — a 

X 



we can write down the exact expressions for A and Ii as 



1 



The two basic quantities that have to be computed in order to solve the dynamic equations for fii and /Z2 are the 
acceptance rate of the Monte Carlo updating: 

(A.6) 



(A.7) 
(A.8) 
(A.9) 

(A.IO) 

(A.ll) 
i: 

(A.12) 
(A.13) 

(A.14) 

(A.15) 
(A.16) 



erfc(a) 



/ 



(ix_ 

erfc(a) +11 — 



2a2 ) px 



The Monte Carlo equations of motion for /ii and /i2 are formally the same found for the fragile glass case 



JQ j dx X W{l3x) p(x|toi, 7712) — (1 + DQ) i dx y{x) W(lix) p(a;|mi, 7712) 



/i2 = / dx X W{j3x) p(x|toi, TO2) + 2/xi / dx y{x) W{l3x) p{x\mi, 11112) 



From (2.23) we know y{x) and we can rewrite it as a function of the above defined a: 

y{x) = 4a2^i (1 - I) 

Using this we get 

JQ - 5^^!Mi±^) _ 4.2,,(i + nQ)Ai^t) 

A^A 



A<2 



K 



hit) + 8a^txiA{t) 
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1. Dynamics in the aging regime: zero temperature 



First of all we solve the equation of motion for /i2 at T = 0, neglecting terms of order /if with respect to those of 
order /i2. For long times is a ^ 1. We can, then, expand Ii{t) for large a, getting 



Equation ( A.16| ) becomes then 



otherwise written as 



hit) 



2ay/7T Aa^ 



(A.17) 



(A.18) 



a ' 



(A.19) 



exp 



At r = 0, the solution in the aging regime, expressed in fi2 — A^/(8q;^), turns out to be: 

A2 1 



log4| + iloglog:% 



\/^F ' 2 ""6 """6 ^ 

Always at zero temperature the leading order of the expansion of the acceptance rate A is , for a ^ 1, 



(A.20) 



(A.21) 



A 



2a0F 



Combining this with ( A.17 ), the Monte Carlo equation of motion ( A. 15 ) takes the form 

e-"' i JQA'^K 



Ml = 



4a2 



-2^ii{l + DQ) (2a^ + l) 



Dividing (A. 23) for (A. IS) we can write down a differential equation for /ii as a function of fj,2' 



(A.22) 



(A.23) 



(A.24) 



where we have neglected terms of order with respect to those of order 1. In the adiabatic approximation, 

obtained by neglecting the left hand side, the solution of (A.24) turns out to be 



Ail 



AJQK 2 
rM2 



A2(l + i^Q)' 

At zero temperature and for long times, one thus have /ii ~ <C A*2. 



(A.25) 



2. Dynamics in the aging regime: T > 

If T is above zero, the leading order of the expansion of A and Ii for large times {a ^ 1) and small temperature 
are: 
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where the terms Ta^ are of 0(1). Indeed it is 



8Ta^ 8T A2 



Ai2 



M2 



(see the definition of a (A.5)) so that 



hm^ = l 



T > 



In this notation (A. 26) and (A. 27) can be rewritten also as 

a" 1 ^ ^ 



A 



h 



P-2 



1 + 



A^k6^2 1 + 



1 + 2 



and the Monte Carlo equations are now 



Ml 



f M2_ 

aV^ 1 + 24^ 

^2 



-4a2(l + DQ)^i 



JQA^K Sfi2 1 

2a2 7I2 1 + 2^ 

M2 



'?A'2 

M2 



aVTT /i2 1^-^ _^ 



M2 



The solution to ( A.33| ) is, to leading order, 



The behaviour of /ii comes out to be 



1 



Ml 



AJQK 



SH2 



M2 



(A.28) 
(A.29) 

(A.30) 
(A.31) 



(A.32) 
(A.33) 



(A.34) 



(A.35) 



Notice that this vanishes when equilibrium is approached, since then Sfi2 0. 

For times even longer than the time scale of the aging regime, the system finally relax, exponentially fast, to 
equilibrium. The equilibrium value of ^2 is known from the statics (see section ^ , as 



T 



M2 = 1> 



(A.36) 



where the explicit expansion of K^{T) in temperature is shown in ( 2.38| ). The asymptotic value of a is, from its 
definition ([A.5|) and taking the first order expansion in T, 



a{T) 



I A2 
8^2 (r) 



(A.37) 
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with 



A'^KD 
8{D + J^) 



(A.38) 



From the equations of motion studied above (look for instance at ( A.19| )) we find for the relaxation time to 
equilibrium 

Teq CX c"' . (A.39) 



Using (A. 37) this is nothing else than the Arrhcnius law 



'oq - exp [ — 



(A.40) 
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